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Abstract 



A multi-channel algebraic scattering theory, to find solutions of coupled- 
channel scattering problems with interactions determined by collective models, 
has been structured to ensure that the Pauli principle is not violated. Positive 
(scattering) and negative (sub-threshold) solutions can be found to predict 
' both the compound nucleus sub-threshold spectrum and all resonances due to 

(f) , coupled-channel effects that occur on a smooth energy varying background. 

f*"^ ■ The role of the Pauli principle is crucial in defining what interaction potentials 

<^^> are required to fit data. The theory also gives an algebraic form for the 

dynamic polarization potential which adds to the ground state interaction to 
' define the optical potential that gives the same elastic scattering cross sections. 

1 Introduction 

o . 

Low-energy cross sections from the collision of nucleons with light-mass nuclei show 
sharp as well as broad resonances upon a smooth, energy-dependent background. 
Those resonances may correlate to states in the discrete spectrum of the target. To 
interpret such scattering data requires use of a complex coupled-channel reaction 
theory . A theory with very important improvements over those used heretofore 
has been developed [Tj. This theory, named a multi-channel algebraic scattering 
(MCAS) theory, finds solution of the coupled Lippmann-Schwinger equations for 
the quantum scattering of two fragments. Though the method can be applied 
to quite general scattering systems, to date studies have been limited to nucleon 
scattering from zero ground state spin targets. 

The prime information sought are scattering (S-) matrices which are easily ex- 
tracted from the T-matrices generated by MCAS. The approach involves using 
matrix algebra on matrices built using Sturmian-state expansions of the relevant 
nucleon-nucleus potential matrix. With this method, all resonances in any energy 
range can be identified, and their centroids, widths, and spin-parities determined 
as can the energies and spin-parities of sub-threshold states of the system. 

It has long been known that collective-model prescriptions of nucleon-nucleus 
scattering violate the Pauli principle. However, the MCAS procedure enables use 
of an orthogonalizing pseudo-potential (OPP) approximation by which unphysical 
effects due to violation of the Pauli principle can be avoided. Doing so is crucial 
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to finding the parameter values of the interaction that simultaneously gives the 
sub-threshold spectrum and the low energy scattering cross sections 

The multichannel formalism and the process by which resonances can be identi- 
fied and located are outlined next. Results of calculations made using a collective- 
model prescription for the potential matrix are then discussed. In that collective 
model, the interaction field is allowed to be deformed from sphericity and that defor- 
mation is taken to second order. In studies of nucleon scattering from 12 C, coupling 
is taken to the ground 0+, l\ (4.4389 MeV), and 0^ (7.6542 MeV) states and we 
consider nucleon energies to 6 MeV. With these three states defining the relevant 
target spectrum, only quadrupole deformation is involved. In a new development, 
we have applied the MCAS method to the p+ 6 He system and have compared results 
with the spectrum of 7 Li known from experiment. 

The theory also allows the construction of the optical potential by appropriately 
subsuming the coupled-channel equations into effective elastic-channel scattering 
equations. The optical potentials thus constructed, as well as allowing for the 
Pauli principle, are very nonlocal and energy dependent. However, the potentials 
guarantee that their use lead to cross section details including resonance effects. 
We present typical results for the n+ 12 C case. 



2 The MCAS theory (in brief) 

The integral equation approach to scattering in momentum space for potential ma- 
trices V^j (p, q) where c = [{£s)j, IJ 71 } where J 71 are conserved quantum numbers 
which will often be omitted for convenience, requires solution of coupled Lippmann- 
Schwinger (LS) equations giving a (partial wave) multichannel T-matrix of the form 



T CC ' (p, q; E) = V cc < (p, q) + p 



open 



/ , i Vcc" (p, x) 72 2 : . Td'd (x, q; E) dx 

.In ki„ - x l + le 



,c"=l 

closed /.oo 9 
/ * , X 



V CC " (P, X) -2 ; 2-T c « c / ( x > T, E) dx 

hin + x z 



(1) 



c"=l 

The open and closed channels contributions have channel wave numbers k c = 
jp(E — e c ) and h c = jfe (e c — E) , for E > e c and E < e c respectively. e c is 
the threshold energy of the target state inherent in channel c and /i is the reduced 
mass. Solutions of Eq. are sought using expansions of chosen model potential 
matrix elements in (finite) sums of energy-independent separable terms, 

N 

Vcc'(p,q) ~ ^Xcnip) Vn 1 Xc'n(q) ■ (2) 
n=l 

In the MCAS method I], the form factors are specified in terms of a set (n) of 
Sturmians $ C)IJ ,(p) via 

\Xc,n)=J2 V cc'\®c,n) , (3) 

c' 

and the Sturmians themselves are solutions of coupled-channel (Schrodinger) equa- 
tions for the chosen matrix of potentials V cc >(p, q) with eigenvalues rj n . 

In coordinate space those potentials are defined using a collective model prescrip- 
tion for the nucleon-target system, with local forms V cc i (r). The Pauli principle can 
be satisfied by using an OPP method in the determination of the Sturmians. Those 
Sturmians then are solutions of coupled-channel equations for the matrix of nonlocal 
potentials 

Vcdr, r') = Vcc>(r)5(r - r') + \A c {r)A c {r') 5 C c- , (4) 
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where A(r) is the radial part of the single particle (bound) state wave function 
in channel c spanning the phase space excluded by the Pauli principle. The OPP 
method holds in the limit A — ► oo, but use of A = 1000 MeV suffices. 

The link between the multichannel T- and S- matrices is S cc ' = 5 cc i—i 
which expands in matrix form with the Sturmians as the basis set to 



-lc + 



j 2fin 



N 

E 

n,n'=l 



kcXcn (k c ) ([t]- Gq] *) xcn' (M \fk^> , (5) 



where now c, c' refer to open channels only. In this representation, rj has matrix 
elements (77) , = i] n S nn i while Gq has elements 



( G o) nn > 



2[i 



E / Xcn(x) 
o=l J ° 



k 2 — x 2 + ie 



— Xcn'(x) d,X 



closed „qq 

E / Xcn{x) 
0=1 70 



(G) 



The bound states of the compound system are defined by the zeros of the matrix 
determinant when the energy is E < and so link to the zeros of { | T7 — Go|} when 
all channels in Eq. © are closed. 

The Sturmian expansions and their use have been detailed recently PP and are 
not repeated herein. Suffice it to note that a finite set can be used to make the 
expansion of finite rank. An essential feature in such studies, if very narrow res- 
onances are to be observed, is a resonance finding scheme. Such is also given in 
detail in Ref. Essentially that requires recasting the elastic scattering S- matrix 
(for each J 7r and with the elastic channel taken as c = 1) as 



Sn = 1 — iir 
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E 

nn' — 1 



k Xin(k)- 



1 r\ 2G0T7 



--Xln'(k). (7) 



Here, the elements of the diagonal (complex) matrix if 2 are — ]==8 nn .>. Then the 

complex-symmetric matrix, t]~^Gq r/ ~i can be diagonalized and the evolution 
of its complex eigenvalues £ r with respect to energy, define resonance attributes. 
Resonant behavior occurs when one of the complex ( r eigenvalues passes close to 
the point (1, 0) in the Gauss plane. At the energy corresponding to that point, the 
elastic channel 5-matrix has a pole. 



3 Results using MCAS for n and p + 12 C 

The MCAS approach has many advantages over the more usual coordinate space 
solution methods. Notably with it one can 1) account for the effects of the Pauli 
principle, even for collective model descriptions of the coupling potentials, 2) specify 
sub-threshold states in the compound nucleus, 3) locate all resonances caused by 
V cc >(r) and find their centroid energies and widths in the energy range chosen, no 
matter how small their widths, and 4) specify the complete T- and S'-matrices off- 
shell enabling prediction of all low energy elastic and inelastic cross sections and 
spin observables. 



3.1 A comparative study 

In this first subsection, we present the results from a comparative study of a coor- 
dinate space program with MCAS (with and without taking into account the Pauli 
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principle). We seek to see if, when one performs as exactly the same evaluation as 
possible, the calculations are equivalent. Also, for a typical low energy problem, we 
ask how the Pauli principle influences the results. 

To make a "fair" comparison, we have used a simplified model for the n+ 12 C 
system. We chose the same three target states to define the coupled channels in 
both the coordinate-space (ECIS97 jHj) and momentum-space (MCAS) evaluations. 
For this we take coupling to be effected by a simple rotational model scheme having 
only a quadrupole deformation, P2 — —0.52, on the real potential 3, 



V(r) = -49.92 /(r) + 



h 



6a-V/(r 



xiv; 

I 



f(r) 



( r-2.885 \ 
o\ 0.63 ) 



(8) 

The potential units are MeV and all lengths are in fm. In the MCAS evalua- 
tion, however, the spin-orbit term is reduced to the Is form. Doing so has little 
effect on elastic scattering [SJ. Initial runs of ECIS97 for the simplified model 

were for (lab.) energies from 0.1 MeV 
to 4.0 MeV in steps of 0.1 MeV. Those 
results indicated that there were three 
resonances near 0.7, 2.1, and 3.3 MeV. 
Thus additional calculations were made 
for more closely spaced energies span- 
ning those three values and the results 
are displayed by the filled circles in Fig. 1. 

Figure 1: The n+ 12 C cross section results 
from using MCAS theory with those found 
using the ECIS approach. Details are given 
E |iib (MeV) in the text. 




Therein they are compared with the results of MCAS calculations, using the same 
simplified model and, essentially, the fixed interaction given in Eq. JHJ, but without 
accounting for the Pauli principle. The solid curve is the MCAS result when the de- 
formation (of V cc i) was taken through second order [TJH]. The three resonances seen 
in the ECIS results are present though their energy centroids are slightly shifted. 
But the background on which those resonances lie agrees very well with the ECIS 
result. The second MCAS result (dashed curve) was obtained by limiting deforma- 
tion to first order and now the agreement with the ECIS result is much better for 
both background and the resonances; centroids and widths. 

We draw two conclusions. First when the simplified model is used in exactly the 
same way in finding solutions of the coupled-channel problem using the coordinate- 
space approach ,3 and the momentum-space approach with MCAS 1\ , the scatter- 
ing cross sections agree very well. The smooth background as well as the specific 
resonances that can be generated with ECIS are found with the MCAS run. Those 
were calculations with deformation taken to first order. The second conclusion from 
comparison of the two MCAS results is that, with deformation of -0.52 in the n+ 12 C 
system, expansion of V cc ' need be taken to second order. 

However, all of the calculations ignored the effects of the Pauli principle [2J- 
But in the MCAS approach, the OPP method can be used to ensure that there is 
no violation of the Pauli principle even when a collective model is used to specify 
the interaction potentials. The OPP method constrains the Sturmians used as an 
expansion set are orthogonal to all states in which the incoming nucleon would 
be trapped into an orbit fully occupied by nucleons in the target. Using such a 
conditioned Sturmian function set to solve the MCAS theory of coupled equations 
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gave the cross section displayed by the solid curve in Fig. 2. There is a consider- 
able change in the resonances observed though the background cross section, being 
dominated by s-wave scattering, is little varied. 
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Figure 2: The n+ C cross section re- 
sults from using MCAS theory with (solid 
curve) and without (dashed curve) using 
the OPP method. 
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Figure 3: The n+ 12 C sub-threshold 
(and resonance) spectra found using 
MCAS with (right) and without (left) 
accounting for the Pauli principle. 



The number of resonances drops from three to two in this energy range when 
the Pauli principle is considered. Further, those resonances inherently have different 
underlying structures, as the | + resonance with centroid near 0.6 MeV is of rank 
1 and not 2 when the Pauli principle is considered. (Here by rank we mean the 
number of the state of given J w in 13 C.) Likewise the \ resonance, in addition to 
being shifted in position by some 2 MeV, has a smaller width and is of rank 2. 

The MCAS method can be used with negative energies and so produce the sub- 
threshold compound nucleus bound state spectrum. With the fixed interaction of 
Eq. JHJ), the spectra (including the resonance states) found with and without using 
the OPP are shown on the right and left of Fig. 3 respectively. If the Pauli principle 
is violated, a large number of spurious states result. 

3.2 Resonance finding 

With MCAS, all resonances (spin-parity, centroid, width) can be found without 
requiring use of an inordinately large number of energy values. The procedure is to 
solve the matrix equations for the eigenvalues ( r (Ei) for a (finite) set of energies. 
The energy centroid of a resonance is that value at which the eigenvalue track 
passes close to the point (1,0). Those eigenvalues are complex and their imaginary 
part relates to the resonance width |Tj. A sample set of the eigenvalues for | + 
resonances in the n+ 12 C system generated using this method are shown in Fig. 4. 
The point values of the largest eigenvalue (labeled (1)) has been connected by a 
long-dashed line to guide the eye to link the energy sequence of the results. The 
actual trajectory has a cusp. That is very evident with curve (1), but similar though 
more slight features are evident in the trajectories (2) and (3). That cusp feature 
of the eigenvalue trajectories links to the opening of the 2+ state at 4.4389 MeV 
in the coupled-channel algebra. Note that these eigenvalues have small imaginary 
parts and so the plot is semi-logarithmic. Thus some details of the trajectories lie 
below the graphed range. The largest eigenvalue clearly evolves well beyond the 
unit circle after crossing at the energy of the lowest ^ compound resonance in this 
cross section. That centroid is located at 3.0 MeV. The imaginary part is consistent 
with the width. Next note 
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that the 2 nd trajectory crosses the 
unit circle at higher energy and with a 
larger imaginary part. This coincides 

o -f- 

with a second, broader, ^ resonance 
at 3.4 MeV. The 3 rd and 4 th Sturmian 
trajectories shown track toward the 
unit circle but have not crossed before 
5 MeV. In this way all resonances of all 
the spin-parities considered are always 
found, no matter how narrow they are. 

Figure 4: Argand plots of | + resonance 
finding eigenvalues for the n+ 12 C system. 
The dashed line is the unit circle. 



3.3 Cross sections and sub-threshold states 

Results of our MCAS study of the neutron plus 12 C system for excitation in 13 C to 
~ 10 MeV, are displayed in Fig. 5. The interaction potential parameters used were 

those tabulated in Ref. |2]. In 
this figure, the elastic scattering 
cross section is compared with 
data with respect to the ground 
state of 12 C as zero energy. 
On the same scale, in the right 
hand parts of the figure, we 
compare the experimental with 
theoretical sub-threshold and 



Figure 5: Spectra of 12 - 13 C and the 
elastic cross section for the n+ 12 C 
system. 

resonance states in 13 C. The labels denote values of 2 J and parity. 

That same elastic cross section and MCAS result are shown in the top panel 
of Fig. 6 with the bottom two panels showing what using the S'-matrices give as 
predictions of polarizations at two scattering angles. The resonance structures most 
evident in the cross section are well reproduced by MCAS and those details are 
confirmed in the polarizations. The spin-parities of the resonances also coincide 
with those of known states in the spectrum of 13 N. Notably the prominent | + 

resonance at 2 MeV, and the two | resonances spanning 2.8 - 4.0 MeV, are found 
with very good widths, narrow and broad respectively. In Fig. 7, proton scattering 
cross sections (from 12 C) is shown at two scattering angles. Found simply by adding 
a Coulomb potential to the n+ 12 C system interactions, the results from MCAS 
calculations (of p- 12 C scattering) compare well with the data. 

3.4 The spectra of 7 Li from the p+ 6 He system 

MCAS results for the p+ 6 He system have been found when the target ( 6 He) struc- 
ture is taken as three states; the + (ground), and two 2 + states at 1.78 and 5.68 
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Figure 6: The elastic cross section for n- 
12 C scattering (top) and polarizations at 
50° (middle) and at 135° (bottom). 



Figure 7: MCAS predictions and proton 
scattering cross-section data at two fixed 
scattering angles. 



MeV respectively. A quadrupole deformation was considered and both 2+ states 
were allowed to couple to the ground with the same deformation and between each 
other through second order. To get the results displayed in Fig. 8, the input in- 
teraction potential required a large diffuseness consistent with the target having an 
extended neutron matter distribution. The results on the far left in Fig. 8, were 
found with the deformation set to zero. They show the sub-structural origins of each 
state in the actual spectrum 0] . Again use of the OPP to account for Pauli blocking 
of occupied states was crucial in finding these results. With zero deformation then 
states reduce to those of a bound proton to one of the three states of 6 He. The 
ground state is a 0p 3 / 2 proton bound to the 6 He ground and 1.8 MeV above that 
there is the quadruplet of states formed by such a proton bound to the first excited 
(2 + ) state. That is also the case for the quadruplet 5.6 MeV above the ground and 
due to coupling to the second excited (2 + ) state. The single state remaining is that 
of a Opi/2 proton bound to the ground. With non-zero deformation, these states 
mix and spread to give the theoretical set that compare quite favorably with those 
known from experiment for 7 Li. Evaluating the n+ 6 He system, one would expect 
more sub-threshold states, as is the case with the mass 13 systems. 

However, one also has to alter the OPP 
since there is a (partial at least) new 

r blocking of the Op-shell. The MCAS 

3 - result for the particle- unstable nucleus 

„ 7 He gives one weakly-bound state. A 

~~ , ' 3 " slight adjustment of the interaction 

. 5" chosen would move that state into the 

_ 7 - continuum without seriously changing 

the 7 Li results. 
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Figure 8: The spectra of 7 Li as determined 
using MCAS for the p+ 6 He system. 



4 MCAS and the optical potential 

Assuming a local form for the elastic channel element of the potential matrix in the 
MCAS approach, and with the operator 



A(E)= V -G^(E) 



n 



(9) 
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involving the free Green function excluding the elastic channel, the optical potential 
for elastic scattering takes the nonlocal form pQ with the addition of 

N 

AU(r,r';E) = £ Xln (r) [A(E)] nn , (E) X i n >{r') , (10) 

n,n' — 1 

to Vii(r). Here Xin( r ) are Bessel transforms of the form factors Xin(k) and AU is 
known as the dynamic polarization potential (DPP). At low energies it is this makes 
the formulated optical potential complex (if the energy allows more than one open 
channel), nonlocal, and energy dependent. The DPP has been calculated using the 

MCAS approach and values of them found 
for 2.73 MeV neutrons on 12 C are plotted in 
Fig. 9 for three partial waves, £ = 0, 1, 2. The 
central depths of each are indicated on the 
right as are the energy steps of the contour 
lines. As the energy is below the first inelastic 
threshold, the DPP is purely real. But the 
coupled-channel effects make that aspect 
of the optical potential extremely nonlocal. 
The DPP is also very angular momentum 
dependent. Of course for energies above 
4.4389 MeV, the DPP becomes complex since 
the threshold for the first inelastic scattering 
has been exceeded. 



Figure 9: The DPP (MeV-fm" 1 ) for three partial 
waves for 2.73 MeV neutrons from 12 C 



The MCAS approach to analyze (low-energy) nucleon-nucleus scattering is built 
from a model structure of the interaction potentials between a nucleon and each 
of the target states taken into consideration. All resonances (narrow and broad) 
in the cross section within the selected (positive) energy range can be found on a 
background. With negative energies, the MCAS method specifies the sub-threshold 
bound states of the compound nucleus. The scheme allows for a resonance finding 
procedure by which all resonances (and sub-threshold bound) states will be found 
within a chosen energy range. Further via use of the OPP in definition of the 
Sturmians to be used, those functions are assured to be orthogonal to any single 
nucleon bound state that is Pauli blocked. That allowance for the influence of the 
Pauli principle is crucial in finding results that concur well with observation. 
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